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We study the influence of phonon scattering on the single-triplet measurement of two-electron spin 
states in a double quantum dot doped with two electrons. We point out that at low temperatures 
phonon-induced relaxation to the ground state suppresses transitions to doubly occupied singlet 
states which are the source of detectable current fluctuations in this measurement scheme. Thus, 
for a relatively strong electron-phonon interaction present in the system, the two configurations 
display the same noise characteristics. In this way, coupling to phonons reduces the distinguishability 
between the singlet and triplet configurations. Under such conditions, the measurement scheme is no 
longer valid even though the temporal characteristics of the measurement process and measurement- 
induced localization are essentially unchanged. 

PACS numbers: 73.21. La, 72.25.Rb, 63.20.kd, 03.67.Lx 
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I. INTRODUCTION 

One of the most promising proposals for solid-state 
qubit implementation is based on the utilization of the 
spin states of electrons confined in quantum dots (QDs). 
The original idea of coding the qubit in the two states 
of a single electron spini still inspires a lot of interest, 
since relatively long spin coherence times have been re- 
ported, and the experimental techniques for the prepa- 
ration, manipulation, and readout of such qubits are be- 
ing rapidly developed^—. The study of two electron spin 
states in double QDs (DQDs) is a natural extension of 
the problem, which serves to examine two-qubit coher- 
ence and inter-qubit interactions— . Furthermore, to fa- 
cilitate electrical control of electron-spin qubits, two-spin 
encoding has been proposed^ - — , which involves spin- 
singlet and spin-triplet configurations serving as the |0) 
and |1) qubit states. This approach proves to be promis- 
ing as well, as is seen, e.g., in the recent demonstration 
of entanglement between two singlet-triplet qubits^. 

The quantum point contaclji 4 - (QPC) measurement of 
charge states in a lateral DQD defined by gate potentials 
in a two dimensional electron gas involves monitoring 
the current flowing through the QPC which depends on 
the occupation of the QDs due to a Coulomb interaction 
between the electrons confined in the QD and electrons 
traveling through the QPC 15 'i2>. This measurement sce- 
nario is a realisation of the so called weak measurement^!, 
where the measured system is only weakly coupled to the 
measuring device. Contrary to the projective measure- 
ment, this measurement is not instantaneous, as both 
the localization of the QD states into the measurement 
basis and acquiring the data needed to distinguish be- 
tween the basis states take time. Apart from the mea- 
surement time, another relevant factor is the attainable 
distinguishability of states, since even after an infinitely 
long measurement time it may not be possible to com- 
pletely distinguish between the measurement basis states. 
On the other hand, a weak measurement is typically 



less destructive to the measured system than an instan- 
taneous projective measurement. Furthermore, such a 
measurement is the only option in many involved quan- 
tum systems which are hard to access experimentally. 
Hence, the QD-QPC measurement setup is commonly 
used experimentally to study QD occupations at very low 
temperatures 1 ^—. As we have previously shown, phonon 
effects do not interfere with the charge measurement in 
any significant way2£, since while they strongly affect the 
coherence times of QD states, phonons do not affect the 
localization times or the distinguishability between the 
measurement basis states in this setup. 

The measurement of spin states of electrons confined 
in QDs is much more complicated and typically involves 
spin-to-charge conversion prior to a QPC measurement of 
the charg o 19 ' 22 ' 24 . A scheme for the direct measurement 
of the spin symmetry (singlet-triplet) of two-electron 
states confined in a DQD was proposed in Ref. [ll|. 
Here, the quality of the measurement relies on QPC cur- 
rent noise being different for the singlet and triplet spin 
symmetries. The disparity of current fluctuations is due 
to the fact that, according to Pauli exclusion principle, 
states with both electrons localized in the same QD are 
allowed in the spin-singlet case, but not for spin-triplet 
states. Hence, the electron charge distribution will fluc- 
tuate during the measurement process due to the QD- 
QPC interaction only if the electrons are in the spin- 
singlet state, leading to enhanced QPC current noise for 
this spin configuration. 

In this paper, we study the interplay of phonon- 
induced effects on two-electron DQD spin states and the 
QPC measurement of these states in the high bias regime. 
In this setup, phonon-assisted interdot tunneling pro- 
cesses at low temperatures lead to relaxation of the spin- 
singlet state to its lowest energy realisation 2 * (with the 
two electrons located in different dots to minimize the re- 
pulsive Coulomb interaction), which should strongly af- 
fect the distinguishability between the spin-singlet and 
spin-triplet states. As we find out, while the additional, 
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phonon-induced decoherence channel obviously increases 
the rate of dephasing of singlet-triplet superpositions, it 
affects neither the time needed for the measured state to 
localize in the singlet or triplet state ( "collapse" of the 
state vector) nor the time needed to perform the mea- 
surement. We show, however, that phonon-assisted tran- 
sitions counteract the enhanced QPC current fluctuations 
in the case of the spin-singlet DQD state by suppressing 
the occupation fluctuations. Thus, the carrier-phonon 
interaction leads to a reduced distinguishability of two- 
electron states. When the strength of the phonon-QD 
interaction is comparable to the strength of the QPC- 
QD coupling, the measurement process is completely sup- 
pressed and an extension of the measurement time cannot 
yield any improvement on the quality of the measure- 
ment. 

The paper is organized as follows. In Sec. [TTl we in- 
troduce the system and define the model to be studied. 
In Sec. IIII[ we derive the quantum master equation in 
Lindblad form for the dynamics of the DQD-QPC sys- 
tem with the electron-phonon interaction included. In 
Sec. IIV1 we introduce the stochastic simulation method 
in the conditional density matrix formalism, which allows 
us to perform simulations of single measurement runs. 
The general results are presented in Sec. [V] while the 
noise characteristics are discussed in Sec. I VII Sec. IVIII 
concludes the paper. 



II. THE SYSTEM AND THE HAMILTONIAN 

We consider two electrons confined in an gate defined 
lateral DQD composed of two identical QDs coupled to a 
QPC, following Ref. 15]. The QPC is located near one of 
the dots (say, right) in such a way that the current flow- 
ing through the QPC is only affected by the occupation 
of this one dot. Since QDs are solid state systems, we in- 
clude the coupling between the confined electrons and vi- 
brations of the surrounding crystal lattice (phonons) , es- 
pecially that the time scales of phonon-related processes 
in such a system are comparable to the times over which 
QPC current traces are observe d 22 i 23 ' 25 . 

The Hamiltonian of the system is given by 
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The first term describes two electrons in the DQD 
structur o 15 i 27 

#dqd = A ^2 {a RrT a La + a\ a a Ra ) +U n i]r n it , 

a=t,i i=R,L 

where A is the amplitude of the tunneling between the 
dots, din-, a) i(y are the annihilation and creation operators 
of an electron in dot i = R, L with spin a = f, |, n ia = 
a\ a a ia gives the number of electrons with spin a in dot 
i, and U is the Coulomb charging energy for adding a 
second electron to a QD. 



I*2> 



E = U 



E = U + J 



WW 

Ci C\ B\ £>2 

\\ \\ 

E = -J 



E = 



\t) 



FIG. 1: Schematic representation of the energy levels of the 
system and the allowed transitions between -Hdqd eigenstates 
induced by the interaction of electrons in the DQD with the 
QPC (d operators, blue arrows) and with the phonon reser- 
voir (Bi operators, red arrows). 



The eigenstates of i?DQD are 
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where \(J(t') — a' La a Ra ,\0} denote singly occupied states, 
and \di) = a^ajjjO), i = L,R, are doubly occupied 
states. The parameters are equal to £ = l/v / 2sin(6'/2) 
and £' = l/V2cos(0/2), where = atan(4A/U). The 
triplet states are degenerate in zero magnetic field, and 
their eigenenergy is chosen as zero. The singlet state 
eigenenergies arc then, respectively, —J,U,U + J, where 
J = l/2(\ / U 2 + 16A 2 — U) is the exchange splitting be- 
tween lowest energy singlet and triplet states (see Fig.rrj. 
The second term describes the QPC lead a 15 ' 16 i 28 ' 29 , 

-Pleads = ^SpC^SpuO-Spa + fr*^2° JD P a Dpo aD P< J > (^) 



p,<7 



where Hu} np is the energy of an electron in lead n = S, D 



npa: a npa 



are the cor- 



(source, drain) and in mode p, a 

responding electron annihilation and creation operators 
with the additional distinction of spin a. The third 
term describes the tunneling of electrons through the 

Q PG 15 i 6 i 29 ) 

-Htun = ^2 (T pq + XpqnR)al pa a Dqa + H.c. (3) 



It consists of two parts: electron tunneling independent 
of the DQD is described by the constants T p9 , while \pq 
quantifies the tunneling dependent on the Coulomb inter- 
action of QPC electrons with electrons in the DQD. This 
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depends on the total number of electrons in the right dot, 
riR = riR-^+nR^. The tunneling constants are assumed to 
be slowly varying over the energy range where tunneling 
is allowe d 15 i 28 , hence we make the assumption T pq s» T 
and Xpq ~ X- We assume that the QPC operates in the 
high bias regime, that is, the chemical potential offset 
between the leads is large enough to induce transitions 
to doubly excited states^ 5 .. 

The last two terms in the Hamiltonian describe the 
energy of the free phonons, 



P h 



(4) 



k,X 



and the interaction between phonons and electrons con- 
fined in the DQD22r— , 



cr,z k,X 



In Eqs (|4| and ([5]), bk,x and bl A are phonon annihila- 
tion and creation operators for a phonon from branch 
A with wave vector k, hujk,x are the corresponding en- 
ergies, F^j' K (k) — F^(k)e ±lkxD ^ 2 are electron-phonon 
coupling constants, and D is the inter-dot distance. 
We include deformation potential and piezoelectric cou- 
plings. The coupling constants for the longitudinal (A = 
I) and transverse (A = tip) acoustic-phonon branches 
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F{k) (6) 



and 
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2p c vw k ,t £q£: 



M tl . t2 (k)F(k), (7) 



respectively, where e denotes the electron charge, p c is 
the crystal density, v is the normalization volume for 
the phonon modes, d is the piezoelectric constant, £o 
is the vacuum permittivity, e s is the static relative di- 
electric constant, and a is the deformation potential con- 
stant. The functions M\ depend on the orientation of 
the phonon wave vector. For the zinc-blende structure 
they are given by 3 ^ 

M\{k) = 2{k x k y (e\ tk ) z + k v k z {ex,k)x + k g k x (ex,k)y], 

where k = k/k and ex,k are unit polarization vectors. 
The form factors .F(fc) depend on wave-function geome- 
try and are given by 



F(k) 
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where tf>(r) is the envelope wave function of an electron 
centered at r = 0. 



III. LINDBLAD MASTER EQUATION 

To describe DQD dynamics averaged over a large 
number of repetitions it is convenient to use the quan- 
tum master equation (QME) approach in Lindblad form 
(Markov approximation). The problem is relatively in- 
volved due to the interaction of our system of interest 
(the DQD) with two reservoirs, a bosonic and a fermionic 
one. Since we assume that these reservoirs are uncorre- 
cted with each other, they can be treated separately, 
yielding unconvoluted terms in the QME2&, 



P(t) 



- [Hdqd , p] 



(8) 
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Here, the first term on the right side of the equation 
describes the free DQD evolution, while the Lindblad 
operators C\ relate to the DQD-QPC interaction, and 
Bi describe phonon-related effects. 

The Lindblad operators (7, may be obtained following 
Ref. [IH in the Born-Markov and rotating wave approx- 
imations (RWA) and assuming independence of the tun- 
neling rates on the initial and final electron state within 
the relevant energy regime. The operators are of the 
formic 



Ci 

c 2 
c 3 



'V-(U + J) . 6. 
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IV + (U + J) . 
I'M z sin-|s Xs 2 | 



(T + ^)I + ivcos-(|siXs 2 | 
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where V — (ps ~ P-d) is the QPC bias, T = ^AngLguT 
is the unconditional tunneling constant related to T vq of 
Eq. and v — \/~§kQl9r is a constant stemming from 
tunneling conditioned on the occupation of the right QD 
(related to Xpq)i where gt is the density of states of the 
i-th lead (i = L,R). C\ and Ci describe inelastic tran- 
sitions which involve energy transfer between the DQD 
and QPC electrons accompanied by transitions between 
the low energy state |sq) and high energy states. The 
quasi-elastic transition between states of similar energy 
\si) and \s2) is represented by the Lindblad operator C3; 
this operator also describes the fully elastic processes cor- 
responding to electrons tunneling through the QPC with- 
out disturbing the DQD state (which are also possible in 
a spin-triplet DQD state). 

To describe the electron-phonon interaction it is con- 
venient to rewrite the appropriate Hamiltonian (Eq. ([5])) 
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in the basis of DQD eigenstates2£, Eqs. (JTJ) , 

ffe-ph = V2[C(\ Sl )( S2 \ + \ S2 )( Sl \) 

-€(l*oX*2|+|flaX*o|)] (h-F R ), 
where the operators are F L / R = Ey^raWfe + 

Following the standard method 38 we obtain the 
phonon Lindblad operators in the Born-Markov and 
RWA approximations, 



B x = 

B 2 = 

B 3 = 
B± = 



/ 72o|s2>(s | 
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y/ 102 e~(u+J)/k B T\ S2){sol 
V721 e-^ T | S iXa 2 |, 



where the transition rates are 7y = 2tt Rij(iVij), with the 
relevant spectral densities of the phonon reservoir defined 
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x [(n k + l)S(w - uJk) + n k S(uj + LJ k )] 
^5>W(fc)| 2 sin 2 ^ 



(9) 



(10) 



k.X 



x [(rife + l)S(w - UJk) + n k S(uj + uj k )] , 

where n k denotes the Bose distribution and Rij{oj) = 
Rji(ui). The energy differences between the states are 
equal to hujQ2 = —h^20 = U + J and Huj 2 i — —hw\ 2 = J. 
Note that phonon related processes involve exactly the 
same pairs of singlet states as transitions related to the 
QPC (see Fig. [1]). In the zero-temperature limit, 720 = 
712 = 0. 

Solving the QME given by Eq. jSJ in the long-time 
limit results in distinct steady states in the spin-singlet 
subspace and the spin-triplet subspace. The triplet 
steady state p^ can be any superposition of the triplet 
states. The singlet steady state is equal to 



P (s) = 

TOO 



Jj [(^++7o2)(^ + 7ai)l«oX*o| 
+ {A 2 o+li2)(A 2 _ +72o)M(si| 
+ (A 2 _ + l2Q )(Al + l21 )\s 2 )(s 2 \] , 



(11) 



where 



N = {A\ +jo2)(Al + 721) + {Al + 7i 2 )(A^ +720) 
+{A 2 _ + 720X^0+721), 



A± = vy/(V ± J ±U)/Hsm(9/2), 
A Q = vy/v/hcos(0/2). 

Hence, the singlet steady state depends explicitly on the 
strengths of the phonon-couplings relative to the QPC 
coupling strengths. 



IV. THE METHOD OF STOCHASTIC 
SIMULATION 

Modern experimental techniques allow one to observe 
single system evolution s 22 ' 23 ' 39 which cannot be repro- 
duced by the QME approach. In our system, this entails 
the situation when, beginning with the same spin singlet- 
triplet superposition state, the DQD system may end up 
in either the triplet state or a singlet-only superposition 
due to the measurement executed by the QPC. Further- 
more, the evolution between initial and final states is 
prone to variation on different runs even if the two states 
are always the same. To describe such a physical pro- 
cess probabilistic elements need to be introduced into the 
evolutioni^^^^-i 2 .. 

To model a single measurement run one introduces a 
conditional density matrix (state) p c (t) that depends on 
the history of counting events (corresponding to electron 
tunnclings through the QPC in our case) and a counting 
process N(t), corresponding to the number of electrons 
that have passed through the QPC. The stochastic equa- 
tion describing the conditional state (in the interaction 
picture) has the form 41 - 

Pc (t+dt) = pc~m)dt[H' cff , Pc u + c ph [ Pc ]dt {i _ dN) 

1 — P\ 



+ EiCiPcC} dN ^ 

"1 



(12) 



where [H e g, p]* — H c sp — pH^, with the non-hermitian 
operator 

ih 



and 



£ P h[p] = 



E U B t B ip + p B l B i) + E B ip B t 



is the Linblad generator accounting for the electron- 
phonon interaction. This dissipative contribution ap- 
pears as a result of averaging over the unobserved phonon 
scattering events and corresponds to the fact that the 
conditional density matrix describes only one subsys- 
tem of the interacting carrier-phonon system. Some of 
the terms proportional to dNdt (that could formally 
be omitted to obtain the equation in its most common 
for m 15 ' 16 ' 29 ' 41 ' 42 ) are kept in Eq. (Ti~2")) for clarity of inter- 
pretation. The increment of the counting process, dN, 
can be zero or one, depending on whether an electron 
tunneling event was observed in the time interval dt. The 
statistics of this increment is defined by the probability of 
a tunneling event (conditional on the history of the mea- 
surement events) Pi = P c [dN = 1] = Tr C\Cip c dt. 

Thus, the first term in Eq. (TT2"j) describes the evolu- 
tion in the absence of electron tunneling events. This 



evolution is described by the deterministic equation 
• /.N, _ Pc - (i/h)[H^,p c }* +£ ph [pc] _ 

Pcwldet — 1 — Pi 

(neglecting terms on the order of (dt) 2 ). This equation 
is satisfied by p c = /9 c /(Tr p c ), where the unnormalized 
conditional density matrix p c evolves according to 

i 

Pc(t)\ det = --[H' cS ,p c }*+£ ph [p c ]. (13) 

Clearly, the trace of p c decreases under the evolution de- 
scribed by Eq. (|T5|) , 

^Trp c (t) = -P 1 (t)Trp G (t). (14) 

The second term in Eq. (fT2"]) , which contributes if a tun- 
neling event has taken place (dN — 1), corresponds to a 
discontinuous change of the system state (a jump). 

The conditional density matrix thus follows the con- 
tinuous evolution described by Eq. (|13[) interrupted by 
jumps corresponding to electron flow through the QPC. 
In order to simulate this piecewise continuous evolution 
we generalize the method of finding the cumulative dis- 
tribution function F(t) for the random jump time t pro- 
posed for stochastic wave function simulations 3 ^ 3 -. We 
note that the survival probability F(t) = 1 — F(t) satis- 
fies the same Eq. (fl"4"|) as Tr p c . Both these quantities are 
equal to 1 at the initial time of the deterministic evolu- 
tion interval. Hence, F(t) = 1 — Tr p c . Based on this, 
the conditional evolution can be simulated by solving the 
deterministic equation for the unnormalized conditional 
density matrix (which, in this case, can in principle be 
done analytically) and finding the next jump time by 
generating a random number according to the known cu- 
mulative distribution function. 



V. RESULTS 

The stochastic method presented in the previous sec- 
tion allows us to model single realisations of the evo- 
lution and subsequently to analyze the current flowing 
through the QPC in a given realisation. In the ab- 
sence of phonons, the noise characteristics of this current 
may serve to distinguish the spin-singlet and spin-triplet 
DQD states—, hence, the interaction between QPC and 
dot electrons can be regarded as a measurement of two- 
electron DQD spin states by the QPC. Since we are 
dealing with a solid state system, one can expect that 
phonon-related effects will disturb this measurement pro- 
cess. 

In the following, parameters corresponding to GaAs 
structures are usecU-^. The DQD energies are U = 1 
meV and J = 0.1 meV, and the QPC bias is taken to be 
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FIG. 2: Exemplary singlet (green, dashed) and triplet (red, 
solid) probability time-evolutions for final triplet (a, b) and 
singlet (c, d) states, without (a, c) and with the phonon in- 
teraction (b, d). 

V = 2 meV. Unless stated otherwise, the QPC tunneling 
parameters are fixed at T = 4- 10~ 2 and v = 9T0~ 4 . The 
material parameters relevant for the calculation of the 
spectral density of the phonon reservoir are^£ cl = 5100 
m/s, c t = 2800 m/s, e s = 13.2, d = 0.16 C/m 2 , a = -8.0 
eV and p c — 5369 kg/m 3 . Two-dimensional Gaussian 
single-electron wave functions were used with 170 nm 
full width at half maximum of the probability density, 
and the distance between the dots was set to D — 250 
nm. This yielded zero-temperature phonon transition 
rates 702 = 1.15 • 10~ 3 ns _1 and 721 = 6.01 • 10~ 8 ns _1 
(720 = 712 = 0). The choice of QPC parameters cor- 
responds to the phonon interaction being roughly 2.5 
times stronger than the QPC interaction, meaning that 

yfito = 2.5Vy/V/H. 

The results presented in this and the following Sec- 
tion are all taken in the zero temperature limit. This is 
because experimental realisations of QPC measurements 
are performed at temperatures that do not exceed 0.1 K, 
leading to extremely low phonon transition rates from 
lower to higher energy states. 

Fig. [2] shows exemplary time evolutions of the proba- 
bility of finding the DQD in a spin-singlet state (green, 
dashed lines) or a spin-triplet state (red, solid) for an 
initial equal superposition state 

l*> = ^(l»> + l*», (15) 

where \t) can be any superposition of the triplet states, 
Eqs l[i"a |) -([i"c ]) . and 

\s) = cos d -\ So ) +sin^| Sl ) = tl> - I It)), (16) 

as in Ref . [l5[ • The top panels show instances where the 
final state is a spin-triplet (the measurement outcome 
was the triplet state) , while the bottom-panel evolutions 




FIG. 3: QPC current as a histogram of tunneling events for triplet (a, b) and singlet (c, d) states, without (a, c) and with the 
phonon interaction (b, d). The time interval used for the histogram is 0.66 /j,s. 



ended up in the spin-singlet state (the measurement out- 
come was the singlet state). The electron-phonon in- 
teraction is included only in the right panels; the phonon 
influence on coherence and localization is studied in more 
detail later on. As seen, regardless of the presence of the 
electron-phonon coupling, the continuous evolution leads 
to DQD localisation in the triplet spin state (see Fig.f2]a, 

b) , while for there to be a measurement of a singlet spin 
state, the occurrence of a quantum jump is required. 

The simulation results for the QPC currents corre- 
sponding to the final singlet and triplet states are shown 
in Fig. [3] Even through the evolutions depicted in Fig. [2] 
show no clear difference between the phonon and no- 
phonon cases, in the phonon-free situation in Fig. [3] 
(left panels) a difference in the magnitude of the cur- 
rent fluctuations (noise) can be seen between the sin- 
glet and triplet case, while no such distinction is evident 
in the current when the phonon influence is included. 
For the realistic choice of material parameters, QPC and 
DQD properties, and for our choice of counting time step, 
the differences are relatively small, but still a period of 
time (after about 110 nanoseconds) when the DQD elec- 
trons occupy higher energy singlet states resulting in in- 
creased current fluctuations can be distinguished (Fig. [3] 

c) . When the phonon coupling is included (Fig. [3] b, 

d) this distinction is diminished (to the level that no 
time period of increased fluctuations can be seen with 
the "bare eye"), so the measurement effect is suppressed. 
The quantitative noise characteristics needed to fully de- 
scribe phonon influence on the distinguishability of the 
spin-singlet and spin-triplet states in the QPC spin mea- 
surement setup are studied in the next section. 

The most obvious effect of coupling to phonons can 
be seen in the dynamics of coherence between singlet 
and triplet states expressed by the normalized ampli- 
tudes of the off-diagonal density matrix elements aver- 
aged over many realisations of single measurement simu- 
lations. The coherence time between the triplet and the 



\s2) singlet states, t S2 = 40.62 fis, is over an order of mag- 
nitude longer than the other two coherence times, which 
are t So = 1.03 fj,s for |s ) and T Sl = 1.20 /us for |si). 
Furthermore, phonon influence (Fig. [5] b) also strongly 
varies depending on the particular coherence in question, 
leaving the triplet \s\) singlet coherence time unchanged, 
slightly influencing the long \s2) coherence, the time of 
which is now t S2 — 40.55 fis, and cutting the coherence 
time of | so) almost by half, leaving t S2 = 0.54 fis. This is 
because different phonon transition rates 7^ govern the 
evolution of different coherences, and the rates signifi- 
cantly differ from one another. The coherence between 
a triplet and the |so) singlet is influenced both by the 
strong transition from \s2) to \sq) and the weaker tran- 
sition from \s\) to Isa), hence phonon influence on the 
coherence times is big. The triplet— |sx) coherence is only 
related to transitions from \sq) to [32) and |s2) to \si) 
which are very weak at low temperatures and do not oc- 
cur at zero temperature, while the \s2) coherence depends 
on the transition rates from \s\) to |s2) and \sq) to \s2), 
of which the first is small and the second is negligible in 
the low temperature range. 

While the dephasing rates are the most common char- 
acteristics of open system dynamics, from the point of 
view of the measurement process the time of localization 
into one of the measurement basis states (the "collapse" 
of the system state) is of more interest. To quantify how 
fast the DQD system reaches the singlet or triplet sub- 
spaces, we introduce a measure of localization analogous 
to the degree of localization used in the description of 
charge measurement^. It determines the timescales on 
which the QPC measurement on the DQD spin states 
is performed. For the singlet-triplet measurement, the 
observable quantifying the localization is (z 2 ), where 

i i 

and the averaging is performed over many simulated mea- 
surement runs. The quantity is equal to zero when the 
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FIG. 4: Localization of the DQD state with singlet post- 
selection (red, solid line), triplet post-selection (green, 
dashed) and with no state post-selection (blue, dotted). 
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FIG. 5: Normalized singlet-triplet coherences without phonon 
influence (a) and with phonon influence (b); |{t|p|so)|: red 
solid lines, |(i|p|si)|: green dashed lines, and |(t|p|s2)|: blue 
dashed lines. Inset: short-time evolution. 



VI. NOISE CHARACTERISTICS 



occupation of the singlet and triplet subspaces is equal, 
and grows with the absolute value of the difference be- 
tween the two occupations, reaching unity when the DQD 
state is either a fully triplet or a fully singlet state. Since 
the singlet subspace is the orthogonal complement of the 
triplet subspace, the degree of localization can be de- 
scribed using the probability of finding the system in 
the triplet state, which gives 2 = 1 — 2(t\p\t). The blue 
dotted line in Fig. 0] shows the localization, (z 2 (t)}, for 
the initial equal superposition state [Eq. (fT3))] for which 
(z 2 (0)) = 0. The red solid line and green dashed lines de- 
pict the evolution of localization for the same set of data, 
on which a post-selection into the set of measurement 
runs that yielded a singlet outcome (red solid) and the 
set that yielded a triplet outcome (green dashed) was per- 
formed. Although the localization occurs on the scale of 
tens of microseconds regardless of the measurement out- 
come, the localization is significantly faster when a sin- 
glet is measured, with localisation time t^ c = 20 /zs, than 

when triplet is measured, with localisation time r£ = 46 
/us. This is because a quantum jump is needed to localize 
in the spin-singlet, and the occurrence of jumps short- 
ens the localization time (as seen in Fig. If no post- 
selection is made then localisation time is t\ oc — 35 [is. 
As expected^, phonons have no effect on the localization 
times as long as they are not monitored (meaning that 
no measurement is performed on the phonon subsystem). 

As can be seen from the above discussion, the phonon- 
induced dissipation destroys coherence but does not 
change the localization time, which depends only on the 
coupling to the measurement device. The major phonon 
effect, from the point of view of the measurement, is the 
reduction of QPC current noise in the DQD singlet state 
due to the phonon induced suppression of transitions to 
the excited singlet states which are responsible for the 
increased singlet current noise. 



Although the destructive phonon effects are already 
seen in the QPC current noise of single evolutions, a gen- 
eral noise characteristic is much more revealing^. Fig. [5] 
shows the two-time maps of the current-current correla- 
tion functioni^^iii, 

G(t,t') = {I(t')I(t))-{I(t')){I(t)), (17) 

where the average (...) is taken over ~ 5000 measurement 
runs (the small-scale fluctuations in the maps are due to 
the finite number of runs). The strongly peaked values 
of the current-current correlation function for t — t' are 
not shown on the maps for the sake of clarity. The ini- 
tial state is chosen to be the equal superposition state of 
Eq. (fl"5|) . In the top panels (Fig. [6] a, b), the correla- 
tion function for post-selected spin-triplet measurement 
outcomes is shown, while in the bottom panels (Fig.[B]c, 
d) the post-selected spin-singlet outcomes are depicted. 
Again, the left panels (Fig. [6] a, c) correspond to the no- 
phonon case and the right panels (Fig.[5]b, d) correspond 
to zero-temperature phonons. When the system reaches 
its steady state, the correlation function no longer de- 
pends on the two times t and t' , but is a function only of 
t = t 1 — t. As can be seen, regardless of the presence of 
phonons, the QPC current noise is Poissonian (meaning 
that there are no correlations present for t ^ t') when the 
DQD is in the triplet state. In the singlet case, current- 
current correlations are present for small time differences 
r, leading to the appearance of the diagonal line (of a 
finite width) on the lower left map. In the presence of 
phonons, these correlations are strongly suppressed due 
to phonon processes that preclude long intervals of occu- 
pation of excited singlet states, and the visibility of the 
diagonal line is much worse. 

The insets on each map of Fig. [5] contain an enlarge- 
ment of the area corresponding to small t and t' , when 
the system has not yet localized in either a singlet or a 
triplet steady state. On the scale of several microseconds 
(consistent with the localization time shown in Fig. 2]), we 
observe small correlations in the upper panels which then 
disappear as the DQD state approaches the purely triplet 
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FIG. 6: Current-current correlation function G(t,t') — G(t = t') of initial equal superposition state for the triplet (a, b) and 
singlet (c, d) measurement outcomes without (a, c) and with the phonon interaction (b, d). Insets: Detail at short times. 



state. On the lower panels we can observe the finite time 
over which the correlations build up, corresponding to a 
decreased visibility of the correlation feature before the 
DQD subsystem localizes in the singlet steady state. 



singlet and triplet steady states, 



To quantify the influence of phonon effects on the 
measurement scheme and its dependence on the ratio 
between the DQD-QPC interaction and the electron- 
phonon coupling we use the Fano factor— ~—. It is a 
widely used noise measure, defined as the zero-frequency 
shot noise power normalized to Poissonian shot noise 
power, F — S(0)/2eI, where / is the mean current. 
The QPC noise power spectrum for steady states is 
given by S(u>) = 2/ dTG(r)e~' lulT , where the current- 
current correlation function, which has been introduced 
in Eq. (|17l) . now only depends on r = t' — t. Following 
Refs [15ll29l | we can calculate the correlation functions for 



G(*/»)(r) = e 2 



-Tr{J2C lP ^Cj} 2 



(18) 



Spectra for the different measurement outcomes can be 
found by substituting the appropriate steady states: the 
singlet state is given by Eq. (ITTj) and = \t)(t\. 

Fig. [7] shows the singlet and triplet Fano factor curves 
as a function of the relative coupling strengths of the 
DQD to the phonon reservoir and to the QPC. For the 
sake of realism it is the tunneling parameters of the QPC 
which are changed, while the electron-phonon interaction 
is kept at a value corresponding to realistic gate defined 
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FIG. 7: Steady state Fano factor for the QPC current as a 
function of the relative strength of the electron-phonon in- 
teraction for singlet (green, dashed line) and triplet states 
(red, solid). Inset: Normalized singlet steady state noise 
power spectra of the QPC detector current without (vio- 
let, dashed line) and with electron-phonon interaction (blue 
solid), a = 1.5. 

QDs (see end of Section |ni] and beginning of Section 
|V|) . The scaling parameter a = ^ = with % = 0.1 
and vq — 2.25 • 1CP 3 , is chosen in such way that a = 1 
corresponds to the situation when the interaction with 
phonons is roughly the same strength as the interaction 
with the QPC; this means that ^702 = v^y/V/h. As can 
be seen, the phonons dominate at large a (small QPC 
current) leading to a suppression of the noise difference 
and breaking of the measurement scheme, while for large 
currents their effect is negligible. Note that the results 
discussed earlier in this paper correspond to the QPC 
interaction strength a — 2.5. Even though the phonon- 
induced suppression of the spin-singlet Fano factor at this 
moderate value of the scaling parameter is small, the ef- 
fects seen in Figs [2] - [5] are already non- negligible. The 
normalized singlet steady state noise power as a function 
of frequency is shown in the inset of Fig. [7J The blue 
solid line corresponds to a = 1.5, while the violet dashed 
line depicts the no-phonon situation at the same QPC 
current values. As can be seen, phonon-interactions lead 
to noise damping for the whole noise power spectrum. 
On the other hand, the width of the noise power spec- 
trum S(ui) does not change when the phonon effects are 
included. This obviously means that the width of the 



corresponding feature in the correlation function G(r) is 
insensitive to phonons. Since this width sets the mini- 
mum time scale over which the measurement has to be 
continued in order to extract the spectal characteristics 
shown in Fig. [7J this leads to the conclusion that also 
this measurement-related time scale is not affected by 
phonon-induced decoherence. 

VII. CONCLUSIONS 

We have studied phonon-influence on the QPC mea- 
surement of two-electron DQD spin states, relying on 
different characteristics of the resulting QPC current 
noise. We have shown that although phonons destroy 
the singlet-triplet coherence, they do not affect the lo- 
calization of the DQD into the measurement basis and 
therefore do not influence this contribution to the mea- 
surement time. Also the time necessary to establish the 
spectral properties of the noise, which is on the order 
of the correlation time, does not change when carrier- 
phonon coupling is included. Nonetheless, phonons do 
disturb the measurement by impeding the distinguisha- 
bility of the spin-singlet and spin-triplet configurations. 
This is due to the phonon-induced suppression of singlet- 
singlet transitions between low and high energy states 
which are responsible for the differences in the noise ob- 
served for the triplet and singlet spin symmetries. 

We have found that the perturbing phonon effects with 
respect to the measurement distinguishability at mod- 
erate relative strengths of the DQD-QPC and electron- 
phonon couplings are already non- negligible. Further- 
more, the electron-phonon interaction can lead to com- 
plete indistinguishability of the two spin configurations, 
if it is strong enough with respect to the DQD-QPC cou- 
pling. This means that the mechanism may render the 
measurement scheme completely useless, while it does 
not impose a lengthening of the measurement times. 
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